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Abstract 

Functional limit theorems are presented for the rescaled occupation 
time fluctuation process of a critical finite variance branching particle 
system in with symmetric a-stable motion starting off from either a 
standard Poisson random field or from the equilibrium distribution for 
intermediate dimensions a < d < 2a. The limit processes are determined 
by sub-fractional and fractional Brownian motions, respectively. 
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1 Introduction 

Consider a system of particles in starting off at time t = 0 from a certain 
distribution (a standard Poisson and equilibrium fields are investigated in this 
paper). They evolve independently, moving according to a symmetric a-stable 
Levy process and undergoing finite variance branching at rate V {V >0). We 
obtain functional limit theorems for the rescaled occupation time fluctuations 
of this system when a < d < 2a. This is an extension of ^ Theorem 2] where 
the starting distribution is a Poisson field and the branching law is critical and 
binary. 

1.1 Branching law 

In the HHHj the law of branching is critical and binary. In this paper an 
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extended model is investigated. The particles branch according to the law given 
by a moment generating function F. F fulfills two requirements: 

1. F'{1) = 1, which means that the law is critical (the expected number of 
particles spawning from one particle is 1 ), 

2. F''(l) < + 00 , which states that the second moment exists. 

(Note here that the branching law in ^ is given by F{s) = 5 (l + s^) and 
obviously fulfills the two requirements.) Although constraints imposed on F 
are not very restrictive and quite natural (so that the class of the branching 
laws satisfying them is broad) still there remain other interesting cases to be 
investigated. One of them is the class of branching laws in the domain of 
attraction of the (1 +/?) —stable law (i.e., the moment generating function is 
F {s) = s + I (1 + s)^'*’^), the case studied in [iC]. A remarkable feature of 
the latter case is that the limit processes are stable ones and not Gaussian as it 
occurs in the finite variance case. 

1.2 Equilibrium distribution 

Another concept naturally related to particle systems is an equilibrium distri¬ 
bution. It has been shown that in certain circumstances the system converges 
to the equilibrium distribution m- It is both an interesting and important 
question whether the theorems shown by Bojdecki et al still hold in the case 
when the equilibrium state is taken as the initial condition. A conjecture in p] 
states that the temporal structure of the limit is given by fractional Brownian 
motion. It is of interest to notice that the limit is different from the one in the 
case of the system starting off from the Poisson field (where temporal struc¬ 
ture is sub-fractional Brownian motion). We study behavior of the system for 
a branching law given by F. But there is still broad area for further studies. 
No attempt has been made to develop more general theory concerning systems 
with a general starting distribution (or a large class of distributions). 

1.3 General concepts and notation 

Let us denote and the empirical processes for the system starting 

off from the Poisson ffeld with Lebesgue intensity measure and the equilibrium 
respectively. For a measurable set A C (A), (A), respectively 

are the numbers of particles of the system in set A at time t. Note that they are 
measure-valued processes but we will consider them as processes with values in 
S' (the space of tempered distributions) because this space has good analytical 
properties. 

The equilibrium distribution is deffned by 

lim 

t —>■ + 130 
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where the limit is understood in weak sense. The Laplace functional of the 
equilibrium distribution is given by 


Eexp{- {Neq,<f)} = exp 


(A,e-^ 


pOC 

l) + v {X,Hiji;s)))ds 
Jo 


( 1 . 1 ) 


where 

j (x,Z) ;= Eexp(-(1.2) 

H{s) = F {s) — s, (f : ^ K+, ip £ £^(]R'^) fl C{TZ'^) and j satisfies the integral 

equation 

j {x, 1) = Tie-^ ix) + V [ Ti_,H {j (•, s)) (x) ds, 

Jo 

This equations can be obtained in the same way as (2.4)]. Note that in ^2] 
function ip is continuous with compact support. We approximate (p G using 
functions ipn with compact support (pn if. Using Lebesgue’s monotone con¬ 
vergence theorem it is easy to obtain the above equations for ip {H is decreasing 
because of the criticality of the branching law). 

For an empirical process Nt the rescaled occupation time fluctuation process 
is defined by 

1 

XT{t) = — {Ns-ENs)ds, t>0, (1.3) 

Jo 

where T > 0 and Ft is a suitable norming. We are interested in the weak 
functional limit of Xt when time is accelerated (i.e., T tends to oo). 

The Qf-stable process starting from x will be denoted by rj^ its semigroup by 
Tt and its infinitesimal operator by A^. The Fourier transform of % is 


ftip{z) = e *1^1 ip{z). 


(1.4) 


For brevity let us denote 


UF (2 - h) 

2 d-i7^d/2Q,p (d/2) h(h — 1)’ 


(1.5) 


where 

h = 3 — d/a (1-6) 

(in this paper we always assume that a < d < 2a so h > 1) and 


M = F” (1). 


(1.7) 


We will now introduce two centered Gaussian processes. One of them is sub- 
fractional Brownian motion with parameter h with the covariance function Ch 


Ch (s,t) =s^ 


(s + t)^ F\s-t\'^ 


( 1 . 8 ) 


and the second one is fractional Brownian motion with parameter h and the 
covariance function Ch 

Ch{s,t) = +t^ - \s (1.9) 
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1.4 Space-time method 


The space-time method is a very convenient technique for investigating the weak 
convergence in the C ([0,r] ,S' (R*^)) space. It was developed by Bojdecki et 
al and can be found in 0. UX = {Xit)\^ [Q is a continuous S' (R'^) -valued 
process we define a random element X of S' (R‘^+^) by 

(^X,<^>)= j\x{t),<^{;t))dt, (1.10) 

where $ G 5 (R‘^+^). In order to prove that Xt converges weakly to X in 
C ([0,t] ,S' (R'^)) it suffices to show that 

(^Xt, (^X,^'j , V$g5(Rd+i) 

and that the family Xt is tight. 

2 Convergence theorems 

We will present two theorems. In the first of them (which is a direct extension 
of m Theorem 2.2]) we study the occupation time fluctuation process for the 
branching system starting off from the Poisson field with Lebesgue intensity 
measure (denoted by A) with the branching law given by a moment generating 
function as described in Section II .11 The result is very similar to the one 
obtained in 0 Theorem 2.2] - namely, the limit process is the same up to 
constants. 

Theorem 2.1. Assume that a < d < 2a and let Xt he the occupation time 
fluctuation process defined by in for the branching system and Ft = 

r(3“a)/2_ Then Xt ^ X in C ([0,t] ,S' (R^^)) as T ^ -|-oo for any t > 
0, where {X (t))(>o is a centered S'-valued, Gaussian process with covariance 
function: 

Cov {{X (s) ,(p) (t ): '0)) = km (A, t) (A, -0) Ch (s, t) , (2.1) 

where v?, '0 G 5 (R'^). 

The second theorem concerns the case where the system starts from the 
equilibrium distribution. As it was mentioned hereinabove the theorem is in¬ 
teresting because the limit has a different time structure from the one in jU 
Theorem 2.2] and Theorem 12.11 

Theorem 2.2. Assume that a < d < 2a and let Xt be the occupation time 
fluctuation process defined by for the branching system V®®, and Ft = 

j,(3-|)/2_ Xt ^ X in C'([0,r],5'(R^)) as T ^ -hoo for any r > 0, 

where {X (t))(>Q is a centered Gaussian process with the covariance function 

Cov {{X (s) ,(p) ,{K (t) , V')) = km (A, ip) (A, V') Ch (s, t), (2.2) 
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where G S {R‘^). 

Remark 2.3. The limit processes above can be represented as follows: 

For Theorem rm 

X = 

and for Theorem 12.21 

where and are respectively sub-fractional and fractional Gaussian pro¬ 
cesses defined in Section o In both cases the limit process X has a trivial 
spatial structure (Lebesgue measure), whereas the time structure is complicated, 
with long range dependence. 

Remark 2.4. The occupation time fluctuation processes of particle systems form 
an area that receives a lot of research attention. We would like to mention some 
other related work. Firstly the case of non-branching systems has been studied 
in 0 Theorem 2.1]. The result is analogous, both to Theorem 12.11 and 12.21 
because the Poisson held is the equilibrium distribution for the system. The 
limit process is essentially the same as in Theorem E2 For the critical d = 2a 
and large dimensions d > 2a, there is no long range dependence and the results 
can be found in . In jH] the fluctuations of the occupation time of the origin 
are studied for a critical binary branching random walks on the d-dimensional 
lattice, d > 3, including also the equilibrium case. The convergence results are 
analogous to those in HE] and in this paper but the proofs are substantially 
different. A similar model with a = 2 was investigated in E] (ie. with particles 
moving according to Brownian motion). 

3 Proofs 

The main idea used in both of the proofs is to study the Laplace functional of 
a process given by the space-time method. The Fourier transform is used for 
this purpose. This is similar to the method in El- In the case of Theorem 12.II 
the proof follows the same principle as E| Theorem 2.2]. The moment gener¬ 
ating function can be represented using Taylor’s expansion and two following 
statements need to be proved. Firstly, one has to check that the method used 
in jH can still be applied. Secondly, it needs to be shown that terms of order 
higher then 2 play no role in the limit. The proof of Theorem I2.2l reniiires more 
work. The Laplace formula contains a function that is a solution of a differential 
equation. This makes the computations more cumbersome. Some expressions 
in this proof had to be examined more carefully than in Theorem 12.1 1 It should 
be noted that Theorem 12.21 covers all branching laws described in Section I lTI 
Now we introduce some notation and facts used further on. 

For a generating function F we define 

G{s) = F{1-s)-1 + s. (3.1) 

The following fact describes basic properties of G which are straightforward 
consequences of the properties of F. 
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Fact 3.1. i. G(0) = F(l) - 1 = 0, 

2. G’ (0) = -F’ (1) + 1 = 0 since F' (1) = 1, 

3. G" (0) = F" (1) < + 00 , 

4-. G (v) = + g {v) where M is defined by il ?. 71) and limj,^o 9 {v) = 0. 

The next simple fact will be useful in proving some inequalities 
Fact 3.2. G (n) > 0 for v G [0,1]. 

Proof. F" (1 — n) > 0 which is an obvious consequence of the fact that all of 
the coefficients in the expansion of F" are non-negative and \ — v G [0,1]. 
G" (v) = F" (1 — n) > 0. We also know that G' (0) = 0 so G' (n) > 0 for 
V G [0,1]. The proof is complete since G (0) = 0 and G is non-decreasing. □ 

The existence of the second moment of the moment generating function F 
implies also that G is comparable with function 

Fact 3.3. We have 

G{v) 
sup — 

«e[o,i] 

Proof. Since both G (v) and are continuous we only have to check that the 
limit of the quotient at n = 0 is finite. This becomes obvious when we recall 
Taylor’s expansion of G (v) from Fact o property 4. □ 

Let us now introduce some notation used throughout the rest of the paper. 
$ will denote a positive function from 5 . [H Lemma in Section 3.2] 

explains why without loss of generality it can be assumed 4) > 0. We denote 

s)= f ^{x,t)dt, 

J S 

{x,s) = jr'i’ ■ 

To make computations less cumbersome we will sometimes assume that 4) is of 
the form 4) (x, t) = tp (x) if (t) for p G S , if G S (K) and hence 

'Srixfi) = (pT{x)xTit), (3.2) 

where pr {x) = ^ip{x), y (t) = J^if{s)ds, xt = x{t)- Notice that p > 
0, X > 0 as 4> > 0. 

Let us introduce now an important function which will appear as a part of the 
Laplace functional of the occupation time fluctuation processes 

(x, r,t) = 1 — Eexp J 4' (•, r -|- s)) , 
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where denotes the empirical measure of the particle system with the initial 
condition Nq = Let us note here that due to the fact that 'f' > 0 we have 
Uii. £ [0,1]. We also write 


nq,{x,rfi)= / Tt-s'^ {■,r+ t — s) {x) ds. 

Jo 

(3.3) 

For simplicity of notation, we write 


VT (x, r, t) = (x, r, t ), 

(3.4) 

riT (x, r, t) = riipj, (x, r, t ), 

(3.5) 

VT {x) = VT (a;,0,T), 

(3.6) 

riT (x) = riT (x, 0, T) 

(3.7) 


when no confusion can arise. 

Now we obtain an integral equation for v which will play a crucial role in 
the next proofs. Note that similar computations can be found also in EH 


Lemma 3.4. satisfies the equation 

vq, (x, r,t) = f Tt-s [4^ (•,r + t - s) (1 - (•,r + f - s, s)) — VG{v^ {x,r + t - s, s))] (x) ds. 

Jo 

(3.8) 


Proof. Firstly let us investigate 


w {x, r, t) = (x, r,t) = Eexp J (TV®, 4'(-, r + s)) ds^ = 1 — (x, r, t ), 

We assume 4' > 0 hence we have w {x, r, t) G [0,1]. By conditioning on the time 
of the first branching we obtain the following equation 

w (x,r,t) =e“^*E y ,r + s)ds'^ 

+ y e“'^®Eexp^-y + u)du^ F {w {r]f ,r + sfi — s)), 

where t > 0, r > 0. 

Using Feynman-Kac formula one can obtain the following equation for w (for 
details see [H (3.13)-(3.17)]) 

■§^w {x, r, t) = (A„ + ^ (x, r)) w (x, rfi) + V [F {w (x, r, t)) - w (x, r, t)], 

w (x, r, 0) = 1. 
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V (x, r, t) = vq, (x, r,t) = l — wq, (x, r, t) so v satisfies the equation 



V (x, r, 0) = 0. 


Its integral version is lfi.81 (note that in m Git) = ^t^). 


v{x,r,t) = I 


Tt-s [i’ {-jr + t - s) {1 - V {■,r + t - s,s)) - VG {v {x,r + t - s,t))] (x) ds. 


□ 


Fact 3.5. 


vq, (x, r, t) < nq, (x, r, t) 


(3.9) 


Proof. This is a direct consequence of the equation the fact that 1 > n > 0 


and Fact o 


□ 


Fact 3.6. For the system the covariance function is given by 


Cov V')) = (A, 


(3.10) 



where G S . 

The proof of the fact follows from a simple computation which can be carried 
on using dni formula (3.14)], therefore we omit it. 

3.1 Proof of theorem 12.11 

3.1.1 Tightness 

The first step required to establish the weak convergence is to prove tightness of 
Xt- By the Mitoma theorem |141 Mitoma 1983] it is sufficient to show tightness 
of the real processes {Xt,4>) for all G <S (K'^). This can be done using a 
criterion ^ Theorem 12.3]. Detailed examination of the proof in ^ reveals 
that only the covariance function of the is needed jH Section 3.1]. One 

can see that the covariance function ifTTn is essentially the same as for the 
binary branching. Hence the proof from ^ still holds for the new family of 
processes. 

3.1.2 The Laplace functional 

The second step uses the space-time method. According to itrm we define Xt 
(from now on r = 1). To establish the convergence we use Laplace functional. 
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By the Poisson initial condition we have (this equation is the same as @1 (3.10)] 


Eexp $^1 = exp ly J (x, s) dsdx^ exp —VTix,0,T)dx^, 

(3.11) 

Now we make similar computations to 0 (3.21)-(3.23)]. By combining II3.1 111 
and lira we obtain: 


Eexp (^Xt, 


= exp < / / (x, s) dsdx 

1 jR'i Jo 


■ exp < - 


= exp 


\ — [ [ 'i’T {x,T — s) {1 — vt ix,T — s, s)) — VG {vt {x,T — s, s)) dsdx 

( Jr'^ Jo 

Iy y 4 't (x, T — s) vt ix,T — s, s) + VG {vt (x, T — s, s)) ds(ia:|' 


dxds. 


The last expression can be rewritten as: 

Eexp{-(lT,$)} =exp{y(/i(T) + /2(T)) + /3(T)}, (3.12) 

where 

Ii{T)=J y ^ ^y {•,T + u - s) (x) dv^ dxds, 

l 2 iT)=[ [ G {vt {x,T - s,s)) - ^ f [ Ty^'i/T{-,T+ u-s){x)du 
Jo Js.'^ ^ \Jo 

(3.13) 

/a (r) = / / ’^T {x,T — s) vt ix,T — s, s) dxds. 

Jo Jr'^ 

To complete the proof we have to compute limits as T —> +oo. We claim 

f f [ ^{x,t)^{y,s)dxdyGhis,t)dsdt, (3.14) 

Jo Jo Jr-^ Jr'^ 

h (T) ^ 0, 

h (T) ^ 0, 

Combining with the above limits we obtain 

lim Eexp (Xt, = exp < f f f f ^ {x,t) ^ {y, s) dxdyCh {s,t) dsdt 

T^+oo L \ / J L 2 Jo Jo Jr'‘ Jr'^ 

(3.15) 

hence the limit process Xt is a Gaussian process with covariance 112.111 . 
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3.1.3 Convergence proofs 


Ii (T) does not depend on F so it can be evaluated in the same way as in ^ 
(3.32)-(3.34)]. 

Let us now deal with I 3 (T). By using Il3.hll we obtain 

I 3 {T)< f ( '^t{x,T — s) f Tu'i'T {■, T — u) dudxds < 

JQ Js.'^ JQ 


C 


T do dR<^ 


ip {x) / Tup (x) dudxds 


Now the rest of the proof goes along the same lines as in 

We will turn to I 2 (T) which is a little more intricate. Combining 113.1311 and 
property 4 from Fact 13.II 


r 

M 

+ u-s){x)dv\ 

+ g{vT{.-.))vT{.-.f 

0 

ft. 

2 

\Jo / 



dxds 


where 


^2 (’^) = [ [ VT {x,T — s, s)^ — ( f Tu"^t {■,T + u — s) (x) dv\ dxds, 

Jo Jr'^ \Jo J 


I2 (T) = [ [ g(vT{x,T - s,s))vt{x,T - s,sf 

Jo Jr-^ 


dxds 


(3.16) 


By inequality (123 we have 

0 < (T) = / / 

Jo Jr'^ 

Combining 123 and (1^ yields 
0 <nT (x, T — s, s) — VT (x, T — s,s) = 


{riT {x,T - s,s)) - {vT {x,T - s,s)y 


[ Ts-u ['I't (•, T -u)vT (', T - u,u) + VG {vT (■, T - u, u))] (a;) du = {*). 
Jo 

We have > 0 for th > 0 which is a direct consequence of the fact that 
T is the semigroup of a Markov process. By Fact 13.31 we have c{F) such that 
F{v) < Hence 


w< 

Jo 

< max(l,c(F)) f %- 

Jo 

< max(l,c(F)) f %- 

Jo 


(•, T -u)vt (•, T -u,u) + c (F) yWT (•, T-U, u)" 


(x) du 


V 


4't (', T -u)vT (•, T -u,u) + —VT (•, T -u,u) 
4't (•, T -u)nT (•, T -u,u) + ^ut (•, T - u, uf 


{x) du 


(x) du. 
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Except of the constant c {F) the last expression does not depend on F. 

Next we consider 

riT {x, T — s, s)+vt (x, T — s,s) < 2nT (x, T — s, s) < 2 f Tg-u^ (•, T — u) [x) du. 

Jo 

The rest of the proof goes along the lines of the proof in 0 inequalities (3.39)- 
(3.42)] and hence we acquire (T) ^ 0. 

Before proving the convergence of (T) we state two facts: 

Fact 3.7. ut (x, T — s, s) ^ 0 in uniformly x € s G [0, T] as T ^ -|-oo. 

Proof. 


Ft 


nr ix,T - s,s) = f Tg-u'^r {-.T - u) 
Jo 

^ lo ^ ( 

/■ 

/ Tu(p [x) du= — 

Jo Jh 


— u) du = 


'T-u 

T 

p{y) 


du < 


x-y\ “ Ft 


0 . 


The last line contains the definition of the potential operator of the semigroup 
Ft which is bounded in respect to x (this can be found in [la Lemma 5.3]). □ 

Fact 3.8. The following convergence holds: 

( f vt {x,T — s,s)'^ ^ c'{'i') as T ^+ 00 . 

Jo Jk'^ 

Proof. One easily checks that 


+ r f VT{x,T-s,sf. 

FI Jo Jm.d 

Hence the result follows from 113.1411 and (F) —> 0 as T —> 0. □ 

It is now easy to prove the convergence of I'J. From Fa,ct l3.1l nronertv 4 we 
know that for given e > 0 we can choose such 5 that Id (^c)] < e. Fact 

K7\ provides us with Tq such that Vt>To (x, T — s,s) < S. Combining this 
with we obtain Vt>To 9 (ft (x, T — s, s)) < e. Hence for T > Tq holds: 


|f^ 2 ^(T)|<e/ f Vt {x,T — s, s) dxds ^ ec {^). 

Jo Jr-^ 

Since e was chosen arbitrary we have convergence: VJ (T) 0 hence also 

h (T) ^ 0 as T ^ -hex). 

Thus we obtained the limits for /i ,/2 and I 3 and the proof of Theorem 12.11 
is completed. 
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3.2 Proof of Theorem 12.21 

3.2.1 Tightness 


We begin by claiming that the family is tight. Close examination 

of 0 Section 3.1] reveals that only the covariance function of the underlying 
system is significant for the proof. By O (3.16)] we know that the covariance 
function of the branching system is of the same form as the covariance function 
of the non-branching system with the Poisson initial condition. From this we 
conclude that Xt is tight. 

3.2.2 Laplace functional for Xt 


We consider Xt defined by irrnt . Using 113 and interchanging the order of 
integration we obtain 


rj' 


UO 


{Nts,'^ i-ys)) ds - { X, J 


To prove the convergence of Xt to X we will use its Laplace functional 

pT ^ 


Eexp{-(i„*.)}= 

{-/ 


'i>T (x, t) dtdx 


Eexpj-/ {Ns,^T{-ys)) ds 


It is easy to check that 
E [ exp i - 


where 


(3.17) 


J (^^s,^'T(-,s))ds| |lVo = = exp{(/i,lnn;T)}, (3.18) 

WT (x) = Eexp |-y (iV“,4'T(-,t)}dt| 


Now we check that 0 < — ln(w 7 ’) is integrable. For T big enough by Fact 13.71 a, nd 
inequity 1S3 we have 0 < ut < c < 1. Hence there exists a constant C such 
that we have — ln(ri;T) = — ln(l — vt) < Cvt < Cut- A trivial verifications 
shows that tit G /1^(R'^) so by HI.HI and II3.18II we obtain 


E exp I — 


{Ns, 5't (•, s)) ds ^ = E E exp 



{Ns,'ifT{-,s))ds 



exp |(A,r(;T - 1) + U J {X, H {Wt {■, s))) ds 
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where Wt satisfies the equation 


Wt {x, 1) = Tiwt {x) + V f Ti^sH {Wt (•, s)) (x) ds 

Jo 

It will be a bit easier to deal with Vr {x,l) = 1 — Wt (x, 1). The equations have 
the form (let us recall that G is defined by lid. 1 1l l 

Eexp|-y (iVs,5'T(-,s))ds| =exp|(A,- vt)+ I^y {X,G {Vr {■, s))) ds^ , 


and 


Vt {x, 1) = Tivt {x)-V [ Ti_sG {Vt (•, s)) {x) ds, 
Jo 


(3.19) 


(3.20) 


Wt is defined by with (p{x) = — lnrt;T(a;) {wt & [0,1] hence ip is positive). 
One can easily see that the definition implies that Wt S [0,1]. Consequently 
Vt S [0,1] which together with Fact 13.21 yields G{Vt) > 0. Hence we obtain an 
inequality 

Vt {x, 1) < Tivt {x) , V3,g]R<i_;>o. (3-21) 

Combining Km and iItT^ we obtain 

Eexp = exp Iy J '^t {x,t) dtdx^ exp J VT{x)dx'^ 

exp-|y y y G {Vt {x,t)) dxdt^ = A{T) ■ B {T), 


where 


A (r) = exp < f f 'S T {x, t) dtdx\ exp ^ — f vt {x) dx\ 


B (T) = exp 


l^y y G{VT{x,t))dxdt} . 


Let us note that A is the same as ll^TI|l in the first proof hence we know that 
its limit is given by 


3.2.3 Limit of B 


To complete the proof the limit limr^+oo B (T) has to be calculated. It suffices 
to consider 

/* + 00 p 




G {Vt {x, t)) dxdt. 


(3.22) 
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Using Fact EH property 4, we split it in the following way 

r [ G{VTi;t))dxdt=^iB,{T) + B2iT) + B3{T)) + B4iT), 
Jo JR"^ ^ 


where 


dxdt. 


BiiT)= / VT{x,tf-{TtVT{x)) 

Jo Jv.‘‘ 

^ + aD p 

B2 {T)= / {Ttvx (x)f - {TtiiT {x))^ dxdt, 

Jo J»^ 

p+aa p 

Bo{T)= / / {Ttur {x)f dxdt, 

Jo Jr‘‘ 

p+co p 

BiiT)= / g{VTix,t))VT{x,tf. 

Jo Jr'^ 

We will prove the following limits (let us recall that we assume II3.2II for sim¬ 
plicity) 

Bi (T) ^ 0, 

B 2 (T) ^ 0, 

BsiT) ^ {X,(pf J j |-Mi - M 2 + (wi + ■“ 2 )^ 1 '0 (ui) V” (^ 2 ) rfui(iu2, 

Bi (T) ^ 0, 

as T ^ -1-00. 

Limit of Bi 

By we obtain 


0 < -Bi{T) = J (^j ^(TtVT{x)f -VT{x,tf dx^ 


dt = 


p + OO p 

/ / {TtVr (x)-V t ix,t)) {TtVT ix)V+T ix,t))dxdt < 

Jo Jr-^ 

Combining this with inequality and equation iOTl we have 

r [ (v [ Tt.t'GiVT{;t')){x)dAi2TtVT{x))dxdt = 

Jo Jr'^ \ Jo / 

Taking into account the form of G (Fact 13.11 property 4) 

Bn (T) + Bi2 (T) , 

where 

Bn{T)= £ £ I ^t-t'VT{-,t' f {x)d£ {2TtVT{x))dxdt, 


14 


Bi 2 {T)= [ f (v [ Tt-t-g {Vt (•, t')) Ft (•, t'f {x) dA {2TtVT (x)) dxdt 

Jo Jv.'^ \ Jo J 

Once again we use inequality 

Sii (T) <VM r [ ( [ Tt-t' iTevT (•))" {x) dA {TtVT {x))dxdt = 

Jo \Jo / 

^ + CX> pt p 

VM / / / Tt-t> {TfVT {x)TtVT {x) dxdt'dt < 

Jo Jo Jw^ 


Applying H^.DII twice 

^+CXD rt 


^+00 pt p 

VM / / / Tt-t'{TfUT {x)TtnT (x) dxdt'dt = 

Jo Jo Jr<^ 


fo JO 

^+00 pt 


P+00 pi p 

MV / / / TfiriT {x)TtmT {x)T 2 t-t'nT {x) dxdt'dt = 

Jo Jo Jk.'^ 


>0 Jo 

We use the Plancherel formula and d 


MV 

~?2d 


I 


+00 pt 


{2tt) Jo Jo Jr^J 

MY c+00 rt 


TtiTiT (zi) TtinT{z2)T2t-t'nT {zi + Z 2 ) dzidz2dt'dt = 


p+oo pt p _ 

/ / / e“* ut (-^i) e~* nT{z2)e ~^‘^^~*ut {zi + Z 2 ) dzidz2dt'dt 

Jo Jo Jr^J 


/ri \2d I / 

(ZTT) Jo Jo 


MV r 

{2n)'^‘^ jR'^'i 

MV f 


nr (zi) nT{z 2 )nT (zi + Z 2 ) 


^ + 00 pt 


to Jo 


^-{2t-t')\z^+Z2r ^^I^MzidZ2 


{2n)‘^‘^ Jr'^J 2 \zi + Z2r (ki|“ + \z2\°‘ + \zi + Z2\°') 
Before proceeding further we will estimate rix 


\nT {z,r,t)\ = 


nr [zi) nT{z2)nT {zi + Z2) dzidz2 = (*) 


Fj' 


f Tt-s'^T (•, r + t - s) ds{z) 
Jo 

[ (p {z) XT F t — s) ds 

Jo 

Bt Jo 


Hence 


|«T {z,r,t)\ < 


C \^{z)\ 


Ft jz 

and this immediately implies (see Il3.7ll i 

C 1 


l_e-*hl“ 


i^T ( 2 :)! < 


Ft ur 




(3.23) 


(3.24) 
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Here, and in what follows, C denotes a generic constant. 
Coming back to (*) and using the last inequality we obtain 


Wl< 


1 


1 


^ ^ 

2 \zi + Z2\ {\zi\ +1^21 + \zi + Z2\ )\zi 




Z2 



1 

X _ g-rlzi+zaT 

L J 

\zi + Z2I 

L J 


Substituting = yi and T^/“z2 = 1/2 yields 


1 


1 


Ct 5 r 

i- I 2 / 1 +2/21“ (|?/ir + | 2 / 2 r + |yi+ 2 / 2 !“) H 
1 


12/2 


l_e-|y2r 

1 

1 _ 

L J 

I 2 / 1 + 2 / 2 I 

L J 


dzidz2 = 


1 - 

dyidy2 < 


where 


(T) ■ B" , 
C’T^ 


(T) = 


F^T^i 


B"i = 


1 


W I2/1+2/21“ (|2/i|“ +12/21“+ l2/i+2/2|“)l2/i| 
1 


1 - 

1 

1 _ p-\yi+V2\°‘ 


I 2/1 + 2 / 2 !“ 

L J 


l_e-l!/d 
dyidy2 


1 2/2 

The integral B'h is finite which will be proved in Fact 14.11 The expression 
B'li (T) can be evaluated 

Bji (T) = T - 5 - = 

and as 1 — ^ < 0 B[^ (T) ^ 0 hence the convergence: Bn (T) ^ 0 is obtained 
too. 

From Fact EH and inequalities (EHl and iHOTll we know Vrix, 2) —> 0 uniformly 
as T ^ 0 and so g (Vr (x, 1)) < e for T sufficiently large hence 

Bi2iT)<e£ (v %_t,VT{-,t’f{x)dt^ {2TtVTix))dxdt<^Bn{T) 

thus Bi 2 (T) 0 and Bi (T) 0 too. 

Limit of B 2 

Let us first estimate expression nr — vt using EH and EH 
nT{x)-VT{x)= [ Tt-u'^’t {•,T - u) (x) du 

Jo 

- [ TT-u['^Ti-,T -u){ 1 -vt{-,T -u,u)) -VG{vTi-,T -u,u))]{x)du 
Jo 
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nT{x)—VT{x)= f Tt-u['^t {-jT — u) vt {'tT — u,u) + VG {vt {'tT — u,u))] (x) du < 
JO 

Applying Fact IlFltl 


Tt-1, 


4't (•, T — u)vt (•, T — u,u) Vc(vt (•, T — u, u))^ (x) du, 


where c is a constant. By inequality 


nT{x)—VT{x)< / "Tt-u — u)nT {-^T — u,u)+ Vc{nT — u,u)) (x) du 

Jo L J 

(3.25) 

We have 


^+oo 

/o 


0 < -B 2 (T) = I (I {TtUT ix)f - {TtVT (x)f dx]dt = 


^ + 00 p 

/ / {Tt {riT {■) - VT {■)) {x)) {% {VT {■) + riT {■)) [x)) dxdt < 

Jo Jrj 

Applying H3.9II and H3.25II 


+ 00 


2 1 'TA Tt 

>0 JrJ Jo 


(•, T — u) ut {■, T ~ u, u) + Vc (jiT (•, T — u, u)) du j {x) 

TtUT (x) dxdt = 

Now we apply the Plancherel formula 


Jo Js. 


Jo 


Is.J' Jo 


'Ft (•, T — u)nT (•, T — u,u) + Vc {ut (•, T — u, u))^ (•) (z) du 

fir (z) dzdt 


Interchanging the order of integration and integrating with respect to t we get 


1 f 1 


(27r) JrJ I^I Jo 


where 


a / Jt—u 


'Ft (•, T — u)nT (•, T — u,u) + VcijiT (•, T — u, u))'^ (•) (z) du 

UT {z)dz = c'{B 21 {T) + B 22 {T)), 


B 2 i{T)= [ < [ 7r_„ [4 't (•,?"-u)nT - U,u)] (z)(iu > UT ( 2 ) dz, 

Jrj fI Jo 


B22{T)= [ -Av / '^T-u Vc{nT{-,T-u,u)f {■) {z)du)nT{z)dz. 

Jrj fI Jo 
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We shall compute limT^+oo ^21 (T) first. We have 


^2 i(T) = e *UT(-,T-M,u)(z)|nT( 2:)d2. 

The inner convolution can be estimated using inequality and simplifica¬ 

tion llfi.2ll 


'i>T{-,T -u)*nT{-,T -u,u){z) = xt{T- u)ipT{-)*nT[-,T-u,u){z) 


< 


XT {T — u) / [z — x) ut {x, T — u,u) dx 
f \!p {z - x) If {x)\ J^dx < 

j-^d |a;| rrp 




In the last inequality we use the fact that f is bounded and is integrable. 
Hence we have inequality 

C 


p2 

Tt 


'I't (',T — u) * riT {• ,T — u, u) (z) < 

Thus B 21 satisfies 

\B2i{T)\<-^[ [ e-^'^-'^^\^\°‘du-nT{z)dz < 

•tR'* 1^1 JQ 

Using inequality and integrating with respect to u 


C 




1 1 


l_e-Ud“ 


1 


J z 


l_e-Ud“ 


dz = 


Substituting = y 

rp3 
-yf ^ 


C 


F^Ti 



1 p-\vr 

1 

1 p-\v\°‘ 

L \yr \yr 


|yr 



(3.26) 


dy < H'l (T) • B" , 


where 




B '21 (T) = C" 


rjrZ 

F^Ti 


f ^ ^ 

1 p-\y\“ 

1 

1 p-i^r 

L \yr \yr 


\yr 



dy 


Is is clear that integral B21 in the last expression is finite since in a neighborhood 
of 0 the integrated expression is proportional to and it is O as 

\y\ -boo (recall that a < d < 2a). Now only B 21 needs to be evaluated 


6 - 3 ( 3-#)-24 -3+^ 

B21 (T) = C"T -3- = C"T^^ 
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Hence it is obvious that B '21 (T) ^ 0 as T —> 0 and so limT^o B 21 (T) = 0. 
Before proceeding to B 22 we will make the following estimation using inequality 


{nr {■,T - u,u)) 

c r 1 

Ft I®I 


{z) = 


/ nr {x, T — u,u) nr {z — x,T — u, u) dx 
JR'I 


< 


1 — e 


l\x\° 


J \z-x\ 


1 — e 


l\z — x\° 


Substitution xv^l°' = y yields 
^ d C 


/ A. 

1 - 

1 [■ 

L \yr 


|zui/“ ~ y\°‘ 


1 — e 




dx < 


, C' 2-^ 
dy < -^u - 


since the integral can be regarded as a convolution of functions so it is 
bounded. This clearly implies 

\B22{T)\<^ [ ^ r e-^^-^^'^'\^-^du-\nT{z)\dz< 

F^ jRd |z| Jo 


C 


„T2 


1 


F^ 


e-(^-“)l"l“dw |nT iz)\dz< 


Using inequality we obtain 

1 1 


C" 


T2-| 

/ x 


Izrizri^ ® ) lzl“ e^l"l)dz 


Fj v/R-i U| Ul 
Substituting zT^^°‘ = y we can rewrite the last expression as 


C 




1 1 


FjTi Vr. \yr i2/r 




The integral is finite (the same proof as for BJ^) and 


T^-i 

FjTi 


= T- 


0-2^ - 3(3- ■^) - 2i 


= T— 


which yields B 22 (T) ^ 0 as T —> +00 

Limit of Bo 

Applying the Plancherel formula to Bo (T) we get 

(ztt) Jo JR'^ 


1 



e = 
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2(27r) 


f T^inTiz)fdz = 

jR'i \z\ 


,d / t\^( f e “)l^l“^r(2)XT(r-M)dM^ dz = 

2(27r)''jR-i |z| \Jo / 


1 1 r 1 


2 (27r)'^ Ft 
Substituting u' = u/T 


f ( f e “I^I“^(z)xt(m)c?w^ dz = 

Jr'^ \z\ \Jo J 


1 f 1 


2 (27r)‘^ -Ft Jr'‘ \^\ \Jo 

1 y2 /•! /•! 


1 


Iff ie-^^“^+“^)l^l“(<?(z))2x(wi)x(«2)dwiciw2dz = 
Jo Jo Jr'* kl 


2 (27r)‘^ 

Let z = [T(mi + M 2)]~“ y 

Z-^—d-T^f f f i'^i+U2)-r^eFy\°‘ (lp([T{ui+U2)]~^y 

2 { 2 Trf Ff Jo Jo JR<i I2/I ^ ^ 

_ ± 

(^ 1 +^ 2 ) x{ui)x{u2)duidu2dy 

Therefore by Lebesgue’s dominated convergence theorem we obtain the limit of 
Bz{T) 


lim B 3 (T) =- ^ 

T^+00 2 (27r) Jo Jo 


(ui+U2)^ “ 1^1 {if{0jfx{ui)x{u2)duidu2dy 

\y\ 


^ [ [ iui+U2)^ “ X (ui) X {U 2 ) duidu2 = 


2'^ar(f)7r2 
Integrating by parts 


0 JO 


K 

W 


(A, (/?)^y J ^-Ui - U 2 + {ui + U2)^'^ tp (ui)'ijj {U 2 ) duidu2 


Limit of B 4 

Firstly, let us notice that 


and hence 


P + OO p 

Bi (T) + B 2 (T) + ^3 (T) = / / Vt (x, t)" , 

Jo JR'^ 

|■+oo p 

/ / Vt ^T^+oc C. 

Jo JR"^ 
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Secondly by Fact I1F7I and inequalities 113/21 II and ll3.DII we know Vr (x) 0 

uniformly as T —> 0. Hence g {Wt (x)) < e for T sufficiently large so 

\B4iT)\<e / VT{x,t)\ 

Jo JR-i 

which clearly implies that B 4 (T) —> 0 as T —> +00. 

Putting the results together 

Combining the previous results we conclude 

^lim B (T) = exp | ^ (A, i/?)^ J J ^—Ui—U2+{ui+U2)^'^tp{ui)ip{u2)duidu2 

And finally by II3.1 511 


lim A(r)B(r) = exp<^ —— (A, / / Ch iui,U2) ^ (ui) tp {U2) duidu2> 

T^+00 I J Jo Jo J 

where Ch is the covariance function of fractional Brownian motion defined by 
oi- This Laplace functional defines a process Xt corresponding to the Gaus¬ 
sian process Xt with the covariance iOl hence Theorem o is proved. 


4 Appendix 

The appendix contains a technical fact used in the main proof. 

Fact 4.1. 


1 


1 


/r- |yi + y2r(|2/ir + |y2r + |2/i + y2r)|yi'“ 

1 


l_e-|yd“ 



1 

1 _ p-\yi+V2\°‘ 

L J 

12/1+2/2! 

L J 


12/2 


Proof. Substituting x = yi + y 2 and z = j/2 we get 


dyidy 2 < -bcx) 


1 


1 


k.. \x\‘^{\x\‘^ + \z\^ + \x-zr) |cr-z|' 


1 _ e-l“-^r 




f —— 

/r- kr i^r 


Let us investigate now 

f 1 


1 


1 

J 1^ 

1 


iQ , I la 


a I iq: 


Rd \x\ + \z\ -I- |a; - z| \x - z 


1 


J X 


1 


1 — 


J ^ 


/Rd |xr + izr + ia;-zr 


1 -e-l" 


J z 




dz < 


f 


1 


, f 1 1 

dz c 1 - -j 


/R-i k! \z\ 


k - z\ 


jR-i k! kl 

L -1 


dz 


dxdz = 

dzdx = (*) 
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The last integral is finite since in the neighborhood of 0 the integrated func¬ 
tion is O and for big Nl is O Going back to (*) we obtain 


(*) < C2 


f —— 

V i^r kr 


1 


-e-M 


< C3, 


by the same reason as above. □ 
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